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■ RESEARCH MEMORANDUM 

LONGITUDINAL FREQUENCY-RESPONSE CHARACTERISTICS OF THE 
DOUGLAS D -558-1 AIRPLANE AS DETERMINED FROM 
EXPERIMENTAL TRANS LENT -RESPONSE HISTORIES ■ 

TO A MACH NUMBER OF 0.90 
By Ellwyn E. Angle and Euclid C. Holleman 


SUMMARY 


During flight tests of the Douglas D-558-I research airplane, 
several transient oscillations of the airplane in response to elevator 
pulses were obtained over a Mach number range of 0-52 to 0-90 at alti- 
tudes between 30,000 and 37,000 feet. ' Through an application of the 
Fourier transform to the experimental input and output functions, the 
longitudinal frequency response of the airplane was determined as a 
function of Mach number . A comparison of the response data estimated 
from wind-tunnel data with the experimental results showed good agree- 
ment. It was found that the maximum response amplitude was a minimum 
at a Mach number of 0.88. At lower Mach numbers (0.52 -to 0 . 66 ) the 
effects of lift coefficient on frequency response were indicated. 


INTRODUCTION 


In the analysis and synthesis of automatic stabilized dynamical 
systems such as the airplane and autopilot combination, there are two 
important elements: ( l) the transfer function of the airplane and 

(2) the transfer function of the autopilot. ■ 


There' is a definite need for airplane frequency- response data at 
high subsonic, transonic, and supersonic speeds to aid in the design of 
automatic-control equipment and to aid in the study of the dynamic 
stability of the transfer function of the airplane. In satisfying this 
need It is more feasible to obtain frequency- response data' from tran- 
sient oscillations in response to pulse disturbances than from the 
response to sinusoidal inputs of various frequency magnitudes . 
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The purpose of this paper is to present the longitudinal frequency- 
response characteristics determined from a Fourier analysis of transient 
oscillations of the Douglas D-558-I research airplane. The data pre- 
sented herein were- obtained over a Mach number range from 0.52 to 0.90 
at altitudes between 30,000 and 37,000 feet-. 
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angle of pitch, degrees or radians 

rate of change of pitch angle, radians per second 

control deflection, degrees 

phase angle between 8 and S e , degrees 

angle of attack, degrees 

stabilizer incidence, degrees 

exciting frequency, radians per second 

undamped natural frequency of airplane, radians per second 
damped natural frequency, radians per second 

damping ratio of airplane tiSl-Sping) 

manipulation variable of Laplace transform 
differential operator per second (d/dt) 
time, seconds 

disturbance-function parameters of transfer function 
Mach number 

normal acceleration, g units 
pressure altitude, feet - 
normal-force coefficient- of- airplane 

rate of change of lift coefficient with elevator deflection, 
per. degree 
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"m p 




j ma 


Jm Da 


'“do 


E e^e 


R 5 e J^5 e 

Subscript: 


■ rate of change of pitching-moment coefficient with 
elevator deflection, per degree 

rate of change of lift with angle of attack, per degree 

rate of change of pitching -moment coefficient with angle 
of attack, per degree 

rate of change of p.itching-moment coefficient with 
angular velocity of angle of attack, per- degree per 
second 

rate of change of pitching-moment coefficient with 
pitching velocity, per degree per second 

real and imaginary parts of output function, respectively 

real and imaginary parts of input function, respectively 


max 


maximum condition 


EQUIPMENT AND PROCEDURE 
Airplane and Instrumentation 


The test airplane, D-558-I, is a high-speed research airplane of 
conventional design with an unswept wing and tail. The physical charac- 
teristics are presented in table I and a sketch of the airplane is 
shown as figure 1. 

All quantities necessary for the analysis of the airplane motion 
were recorded by standard NACA recording instruments and were synchro- 
nized by a common timer. The quantities recorded were: normal accelera- 

tion, airspeed, altitude, pitching velocity, elevator position, and 
stabilizer position. Accuracies of the recorded quantities are indi- 
cated by the following instrument accuracies: 


Pitching velocity, 6, radians per second . . . ±0.005 

Elevator position, S e , degrees ^ ±0.1 

Mach number,. M ....... • ±0.01 

Normal acceleration, n, g units . . . ±0.05 
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Method of Obtaining Data 

During flights of the D-558-I airplane , several stick-fixed 
longitudinal responses to elevator pulses were obtained. The responses 
were recorded during stabilized 1 g flight over a Mach number range 
of 0-52 to 0.90 at— altitudes between 30,000 and 37,000 feet. The input 
function consisted of a pulse of the elevator of about 2° in magnitude 
ahd a duration of from 0.5 to 1.0 second. 

The pulse resulted in an initial airplane oscillation of approxi- 
mately ±l/2 g to ±1 g and a pitching velocity of ±0.1 radian per second. 

A restricting device was used on the elevator control to insure the 
return of the control to approximately its original position following 
the disturbance and to aid in maintaining fixed-elevator condition 
during the subsidence of the oscillation. All other controls, ailerons, 
rudder, ■ and stabilizer, were fixed during the maneuvers . The stabilizer 
was fixed at -2.00°, airplane nose down, throughout _ the tests. The 
elevator position was recorded at. four stations along the span external 
of the fuselage. In order to compensate for elevator twist, which 
amounted to about 0.4°, an average of the four- control positions was 
used as the input function. 


METHOD OF ANALYSIS 


The analysis of the transient- response flight data Is based on the 
assumptions that the equations of motion are linear differential equa- 
tions with constant' coefficients and that-forward velocity and altitude 
are constant. The equations of- motion are as follows: 


c Ls 5 e + c L a a = 2 t D ( a " e ) 


( 1 ) 


C m$- 5 e + ' C m a a + C m Da 1)0 + C mpe 1,6 - 2V 


( 2 ) 


■ where 


t = 


1 m . 

pSV 


K 2 - kly 
- 2 

me 
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m 

S 

V 

p 

Iv 


mass of airplane, slugs 
wing area, square feet 
forward velocity, feet per second 
air density, slugs per cubic foot 

r 

moment of inertia in pitch, slug-f eet £ 
mean aerodynamic chord, feet 


Time histories of pitching velocity and elevator position provide - 
the data from which the transfer functions of the system are "determined 
by means of the Fourier transform. The use of this transform is briefly 
described herein but a more complete explanation may be found in 
references 1 to 3- 


The Fourier transform, as applied in harmonic analysis, is 

-tot 


F( to) = J 


F(t)( 


dt 


(3) 


By transferring the experimental functions , input and output, so that 
the initial and final conditions are equal tqpzero, the transform may 
be simplified to 

F(to) = [ Fftje' 1 ®* dt (4) 

J 0 

where the quantity F(t) reaches a steady- state condition at .time t = T 
identical to the initial condition at the time t = 0. If,' however, 

F(t) assumes some steady-state condition other than that equivalent to 
the initial condition, the transform is modified to 

t 

F( to) = f F(t)e" tot dt + F F(t)e _:kDt dt ■ (5) 

Jo . Jt 

for which the second term on the right-hand side becomes finite as 



( 6 ) 
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Evaluation of the integral equation (5) yields a set of real and 
imaginary parts each for the input' and outputr-as a function of consecu- 
tive integral multiples of the fundamental frequency cuq- The integrals 
in the real and imaginary parts of equation (5) are obtained by numerical 
integration (Simpson's One-Third Rule) of 1 the respective parts as 


R^ = r 0(t) cos cut dt - ^ sin off 


pT 

= - / S(t) 

W A 


sin cut dt - — ^ cos ctT 


cu 


(7) 


R 6 e 


r T 

= / 5 e (t) cos cut dt - 

J o 


5_(T) 


cu 


sin off 


x 5 ■ = - / ~e 
e J 0 


, X Sg(T) 

S e (tr) sin cut dt - — cos afC 


cu 


including the corrections f-pr a finite steady-state condition. 

It is possible to combine these quantities directly into a 
frequency-response expression as 


Q R0 +-ile 

6 e = R 5e + H 6e 




(8) 


The results determined by equations (7) and (8) are exact within the 
limitations of numerical integration accuracy. 


The frequency response can be determined in- another way from 
transient-response data by assuming a form for the transfer function as 


Cq + c^s 


( s2 + + 2 Ia% 


( 9 ) 


derived from the assumed two -degree -of-freedom system described by 



NACA EM L51K28 


7 


equations (l) and (2) . These results, however, can only be classed' as 
faired data. Application of the incomplete Fourier transform to equa- 
tion (9) yields: 


(icu)^ + 2|^a>g^iaj + J 0(t)e - ' im " ,: ' dt + e*^ 


j ^ 0 


[d0(T) + 21^^ (T) 


i£u6(T)J - jneCo) + 2| A 05 RA e(0) + iux9(0)J = (Co 4 

Cjico) [ SgCtJe^ dt + C 1 5 e (T)e' cjjr _ 0^(0) 
w r\ "" 


( 10 ) 


from which a real- and an imaginary equation can be set up. By using the 
method of least squares and a matrix solution, the transfer.- function 
coefficients Cq, Cj_, 2! a o> r ^, and are determined from the solution 

■of the simultaneous equations derived from the number of frequencies 
desired to give a frequency spectrum. (Least- squares method is described 
in reference 4. ) The transfer function for the system is then calculated 
for the desired exciting frequencies by substituting the transfer- 
function coefficients back into equation (9) • 


RESULTS AND DISCUSSION 


Presented as figures 2(a) to 2(d) are four representative time 
histories of the quantities and the type of transient oscillations used 
in the determination of the transfer functions. The numerical trans- 
formation of time histories such as these from the time plane to the 
frequency plane yields the frequency- response data shown as figures 3(a) 
to 3( <0 which show response amplitude and phase angle plotted against 
frequency. 

The points in the figures represent the experimental frequency 
response determined by equations (7) and (8) . The same real and 
imaginary parts, equation (7) without end correction, are used in con- 
junction with equation (10) and the least-squares method of curve fitting 
to give the variations represented by the continuous lines and presented 
as faired data. These results show the variation of the frequency- 
spectrum characteristics at constant Mach numbers over a Mach number 
range of 0*52 to 0.90. It can be seen that the maximum response ampli- 
tude and the natural frequency take on a definite variation with, Mach 
number. This variation is demonstrated in figure 4 in which these 
quantities are plotted against Mach number. The maximum amplitude 
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decreases from 9-3 at a Mach number of 0.52 to a minimum of 6.4 at a 
Mach number of 0 . 6 ^ } increases to 9-9 at a Mach number of 0.75> then 
decreases to a minimum value of 2.7 at~a Mach number of- 0.88 after which 
the curve again increases to a value of 7*1 at— a Mach number of 0.9- 
Over the same Mach number range the natural frequency decreases from a 
value of 4.00 at a Mach number of 0.52 to 1.7 at a Mach number of 0.75 
and then increases to a value of 5-2 at a Mach number of 0-90. 

Between Mach numbers of 0.60 and 0.90 the maximum response amplitude 
and natural frequency derived from wind-tunnel data (reference 5) acre 
compared with the experimental results determined from Fourier analysis. 
The agreement is considered good from Mach number 0.69 to Mach number 0.90 
The decrease in the natural frequency might possibly be attributed to a 
lift-coefficient effect since the lift coefficient decreases from 0.8 
at a Mach number of 0.52 to 0.4 at a Mach number of 0-75* 


CONCLUDING REMARKS 


An application of the Fourier transform has been used to obtain 
the frequency-response characteristics of- the Douglas D-558-I airplane 
from experimental, transient-oscillations in response to a pulse elevator 
input over a Mach number range of 0.52 to 0.90 at altitudes between 
30,000 and -37* 000 feet. The experimental frequency- response data were 
obtained by numerically applying the Fourier transform to the transient 
response obtained. A comparison of the response data estimated from 
wind-tunnel data with the experimental results showed good agreement. 

It was found that the maximum response amplitude had a minimum value at 
a Mach .number of 0.88. At the lower Mach numbers (0.52 to 0.66) the 
effects of lift coefficient on the frequency response were indicated. 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I . 

PHYSICAL CHARACTERISTICS OF DOUGLAS D- 558-1 AIRPLANE 

Airplane : 

■Mass, slugs . . . . 315 

Moment of-inertia In pitch, slug-ft 2 l6, 750 

Fuselage: 

Length, ft 35-04- 

Wing: 

Span, ft . . 25 >0 

Area, sq ft 150. 7 

Taper ratio ....... 0.54- 

Aspect ratio . . . . •• 4.17 

Mean aerodynamic chord, ft 6.21 

Horizontal tail: ' 

Span, ft 12.25 

Area, sq ft 35-98 

Taper ratio 0-55 

Aspect ratio 4.17 

Tail length, from 0.25 chord line of . wing to hinge line, ft . . 16.34- 





the Douglas D-558-I airplane 
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T/ me j t , sec 


(a) M = 0.52; h p = 37,400 feet- 

Figure 2.- Time histories of. typical airplane response to an elevator 

pulse, i-fc = -2.00°. 
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(t>) M = 0.55;. hp =’35,800 feet;. 
Figure 2.- Continued. 
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7~//??&, tj sec 


(a) M = 0 . 90 ; hp = 30,000 feet. 


Figure 2.- Concluded. 







ratio > ejd t 



O Expen/wenta/ dj8 e 
a Experrmerta/ ft 
~ Fa/red a/5 e 
Faired <p 


Frequency >(*}> rad/ a ec 


(a) M = 0.52; h p = 37,^0 feet. 

Figure 3*- Frequency response of the Douglas D-558-I airplane as deter 
mined from pitching velocity response to an elevator pulse. 
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Frequency > CO , rad/ ■sec 


(c) M = 0.55; hp = 36,400 feet. 
Figure 3>- Continued. 


Phase angle t (f> , deg lead 









4 6 8/0 

Frequency > CO; rad/ sec 

(e) M = 0.65; hp = •40,200 feet. 


Figure 3-- Continued. 
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(h) M = 0:75; bp = 38,000 feet. 
Figure 3-- Continued. 


Phase angle , <f> > deg lead 
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Frequency t CO, radjsec 

(i) M *= 0.82; hp = 3^,000 feet. 
■Figure 3<- Continued. 


Phase ang/e, cp, deg Jead 





/^repuencj/j 6J } rad/vec 


i ' M = 0.8 6 ; tip = 33,000 feet. 


Figure 3-- Continued. 
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frequency ,C0 , rad /sec 



(Z) M = 0.90; h p = 30,000 feet. 


■ Figure 3*~ Concluded. 
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2 


40 



.60 SO 

Mach number, M 


t4 


IZ 


fO 


6 


ro 

oo 


/.00 


Figure 4. - Comparison of wind tunnel and experimental maximum response 
amplitude and damped natural frequency as a function of Mach number. 
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